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Abstract 

The shape of Bose-Einstein (or HBT) correlation functions is deter- 
mined for the case when particles are emitted from a stable source, ob- 
tained after convolutions of large number of elementary random processes. 
The two-particle correlation function is shown to have a stretched expo- 
nential shape, characterized by the Levy index of stability < a < 2 
and the scale parameter R. The normal, Gaussian shape corresponds to 
a particular case, when a = 2 is selected. The asymmetry parameter of 
the stable source, (3 is shown to be proportional to the angle, measured 
by the normalized three-particle cumulant correlations. 

In high energy nuclear and particle physics, the space-time structure of par- 
ticle emitting sources is often investigated with the help of the two-particle 
Bose-Einstein correlation functions. In heavy ion physics, these correlations are 
frequently called as HBT correlations to honor the astronomers R. Hanbury 
Brown and R.Q. Twiss, who invented a similar method pQ in radio astronomy 
to measure the angular diameter of main sequence stars. 

The two-particle correlation function C2(ki,k2) is defined as 

( JV 2 (ki,k 2 ) m 

where kj is the momentum of particle i — 1,2, and A^ 1 (k 1 ) is the single particle 
invariant momentum distribution (IMD), while ^(kijka) is the two-particle 
invariant momentum distribution. 

In this manuscript we highlight some of the results of ref. 2 , where we have 
investigated in great detail the Bose-Einstein or HBT correlation functions under 
the following three 

experimental conditions: 

i) The correlation function tends to a constant for large values of the relative 
momentum q — k% — tc2- 

ii) Near \q\ = 0, the correlation function deviates from its asymptotic, large 
\q\ value in a certain domain of its argument. 
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Hi) The two-particle correlation function is related to a Fourier transformed 
space-time distribution of the source. 

Condition Hi) is satisfied if the propagation of identical boson pairs from a 
chaotic (thermalized) source to the detector can be described by a plane wave 
approximation; this is possible if Coulomb and strong final state interactions as 
well as additional short range correlations e.g. caused by resonance decays can 
be corrected for or are negligible. 

For clarity, let us consider first a one dimensional, factorized toy model. Let 
x and k stand for the coordinate and momentum variables, respectively. The 
model is defined by the emission function 

S(x,k) = f(x)g(k), (2) 

and the normalizations are 



dxf(x) = l, Jdkg(k) = (n), (3) 

where (n) stands for the mean multiplicity. The single-particle spectrum is 

Nx(k) = [ dxS(x,k)=g(k). (4) 



If Hi) is valid, the Bose-Einstein symmetrized two-particle wave-function is 

i/>ki,ks(xi,Z2) = —j= [exp(jfcixi + ik 2 x 2 ) + exp(ikix 2 + ik^xx)] . (5) 
V 2 

In the Yano-Koonin formalism [3], the two-particle momentum distribution of 
chaotic sources is given as 

N 2 {k ll k 2 ) = J dxidx 2 S(x 1 ,k 1 ) S(x 2 ,k 2 ) \ip kl ,k 2 {xi,x 2 )\ 2 . (6) 

Let us introduce as auxiliary quantities the Fourier transformed source density 
distribution and the relative momentum as 

/(?) = J dx cxp(iqx) f(x), q = k\ - k 2 . (7) 

The two-particle Bose-Einstein correlation function is obtained as 

C 2 (k u k 2 ) = l + \f(q) | 2 , (8) 

that measures the absolute value squared Fourier transformed coordinate-space 
distribution function of the particle emitting source. 

In physics, as well as in the theory of probability, the probability distribu- 
tion of a sum of a large number of random variables is one of the important 
problems, and such distributions are frequently realized in Nature. Limit dis- 
tributions characterize the probability distributions of random processes in the 
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limiting case when number of the elementary independent random subprocesses 
tends to infinity. In case of high energy nuclear and particle physics, for ex- 
ample, the position of emission of an observable particle is obtained as a sum 
of a large number of position shifts due to various parton-parton scatterings, 
hadronization, rescattering of hadrons, and decay of hadronic resonances: 

X = ^2 x n (9) 

n 

Hence the distribution of the sum x is obtained as an n-fold convolution, 

f(x) = J dx 1 ...dx n fi(xi)...f n (x n )5(x-xi-X2....-x n ) (10) 

Various forms of the Central Limit Theorem state, that under some conditions, 
the distribution of the sum of large number of random variables converges to 
a limit distribution. In case of "normal" elementary processes, the limit distri- 
bution of their sum is the Gaussian distribution. This is one of the frequently 
encountered cases of limit distributions. 

Stable distributions are precisely those limit distributions that can occur in 
Generalized Central Limit theorems. Their study was begun by the mathemati- 
cian Paul Levy in the 1920's. A recent book by Zolotarev and Uchaikin 4 
contains over 200 pages of applications of stable distributions in probabilistic 
models, correlated systems and fractals, anomalous diffusion and chaos, physics, 
radiophysics, astrophysics, stochastic algorithms, financial applications, biology 
and geology. Stable distributions provide solutions to certain ordinary and frac- 
tional differential equations. The breadth of their applications suggests that 
they can be considered as a class of special functions @J|51|I)]. 

The Fourier transformed density distribution is usually called the character- 
istic function in mathematical statistics. The stable distributions are frequently 
given in terms of their characteristic functions. The reason for this is that the 
Fourier transform of a convolution is a product of the Fourier-transforms, 

n 

/(ff)=n&«) ( u ) 

i=l 

and limit distributions appear when the convolution of one more elementary 
process does not change the shape of the limit distribution, but it results only 
in a modification of the parameters of the limit distribution. Hence the stable 
distributions have simple characteristic functions. However, the explicit for- 
mulas describing the Levy stable source density distributions are known only 
in some special cases. As of now, this is not an essential limitation as public 
domain numerical packages exist that can be utilized to calculate these source 
densities for any set of parameters [7] . 

Results of mathematical statistics yield a simple form for the characteristic 
function of univariate and symmetric stable distributions, 

f(q)=exp(iqS-\jq\ a ), (12) 



3 



where the support of the density function f{x) is (—00, 00). Deep mathematical 
results imply that the index of stability, a, satisfies the inequality < a < 
2. This parameter determines, for large modulus of the coordinates, the Levy 
distributions. Levy laws with index of stability a tend to power-laws and the 
exponent of the decay of these distributions is given by 1 + a, f(x) — > 
for |x| — > 00. 

Although the proof that the shape given in eq. l(T^|) is unique and that 
< a < 2 is rather complicated, it is easy to show that these Levy distributions 
are indeed stable under convolutions, 

n 

Mq) = exp (iqSi - 17,0 , II = ex P ^ 5 ~ W) > ( 13 ) 

1=1 

n n 

r = £7?, * = (14) 

i=l i=l 

so after appropriate shifting and rescaling the stable distributions remain in- 
variant. Observe that eq. I|14l) generalizes the well known quadratic addition 
rule of variances of convoluted Gaussian distributions to the case of stable dis- 
tributions. 

In the following, let us adopt the notation of Nolan Our choice corre- 
sponds to the S(a, = 0, 7 = R/2~ , 6 = xq; 1) convention. In order to simplify 
the results, and to present results that are similar to the ones used in data fitting 
in high energy and nuclear physics, we have re-scaled the scale parameter 7 of 
the Levy distributions and introduced a physical notation as follows: 

R = 2*7, x = S (if /3 = 0). (15) 

In the chosen S(a, (3, 7, 6; 1) notational system, for symmetric ((3 — 0) stable 
distributions, the parameter S coincides with xq, the location parameter of the 
distribution, characterizing the position of particle production. (This parameter 
is irrelevant in Bose-Einstein correlation studies, we shall see that it cancels from 
both the two- and the three-particle correlation functions.) In this notation, the 
one dimensional symmetric Levy-stable distribution yields the following, simple 
form of the two-particle BEC: 

C(q) = l + CX P (-\qR\ a ). (16) 

This form that has an additional parameter, the index of stability a, as com- 
pared to the usual Gaussian (or exponential) distribution, where the value of a 
is fixed to 2 (or 1). This result can be generalized straightforwardly to multidi- 
mensional expanding, core-halo type of systems. Two examples are given here, 
for more throughout discussion see ref. 0- 

Example a). For collisions with non-relativistic energy, and a small duration 
of particle emission, symmetric Levy distributions yield the following Bose- 



4 



Einstein correlation function [2]: 



C 2 (fci, fc 2 ) = 1 + Aexp 



(17) 



As usual, core-halo corrections |Sj introduce the intercept parameter A. The 
three-dimensional expansion of the core results in a multivariate decomposition 
of q. All the fit parameters may depend on the mean momentum, (A, R 2 j, a) — 
(A(K),i??.(K),a(K)), where K = 0.5(ki +k 2 ). 

Example b). For collisions with very high energy, the particle emission pro- 
cess is a highly relativistic phenomena. In this case, the invariance of the emis- 
sion function can be reflected if the longitudinally boost-invariant proper-time 
variable r = y/t 2 — r 2 . is utilized, and the space-time rapidity i] = 0.5 log [(t + 
r z)/{t — r z )\ is also introduced as a hyperbolic, boost-additive coordinate. 

In a factorized form the Buda-Lund (BL) parameterization assumes the fol- 
lowing structure for the emission function [Sj: 

S(x,k) = H4T)G*(r))h(r x ,r y ), (18) 

where the subscript * denotes an implicit momentum dependence. The effective 
proper-time and space-time rapidity distributions iZ* (t) and G* {rf) are assumed 
to have a uni-variate Levy distributions with indexes of stability a = and a.\\, 
while i*(r x , r y ) may have a bivariate Levy distribution with index a±. We use 
the symbolic notation of the invariant temporal and the parallel relative 
momentum differences, Q = and Qm, being conjugated variables to the space- 
time variables (r, rf). In these variables, the correlation function is 

C{kxM) = l + Aexp(-|.R = Q = r= -|i?||Q||| Q ii -\R±Q ± \ a ^), (19) 
Q= = m t ,\ cosh(yi - rj) - m 4:2 cosh(y 2 - rf) (20) 
Q\\ = m t> i sinh(yi - rf) - m t ,2 sinh(y 2 - rf) (21) 



Ql + Ql- (22) 



In these equations, m t i — J m 2 + kf is the transverse mass, yi — 0.5 ln[(£Ij + 
k z i )/(Ei — k z i )] is the rapidity of particle i and the fit parameter rj stands for 
the space-time rapidity of the point of maximum emissivity for particles with a 
given fixed four-momentum ki jlUI 2] . The three different indexes of stability 
satisfy the usual inequality < a.i < 2 for all ; = (=, ||, ||). All the five fitted 
scale parameters, A, R^, R,± and rj, as well as the three Levy indexes a = , 
a\\ and aj_ may depend on the value of the mean momentum K. 

Finally, let us consider the case of three-particle Bose-Einstein correlations. 
If the particle emission is completely chaotic and the plane- wave approximation 
can be warranted, this reads as 

C 3 (l,2,3) = 1 + |/(1,2)| 2 + |/(2,3)| 2 + |/(3,1)| 2 + 

+ 2Hf(l, 2)/(2,3)/(3,l). (23) 



5 



where the symbolic notation f(i,j) = f(ki — kj) = /(%•) has been introduced 
to simplify the equation. The three-particle cumulant correlation function cor- 
responds to the last term, 



K3(l,2,3) = 2ft/(l,2)/(2,3)/(3,l), 



(24) 



where the two-particle cumulant correlation function is defined as 



K 3 (l,2) = |/(l,2)| a . 



(25) 



Hence, the normalized and symmetrized ratio 



w(l,2, 



3) 



«s(l,2,3) 



(26) 



2^(1, 2)k 2 (2,3)k 2 (3,1) 



turns out to be a simple function of (3, the asymmetry parameter: 



Ml, 2,3) 




fiTtan(— )Q2)\<ki\ a *g*(<Hi)]\ fOT (28) 



(for the special case of a = 1, see again ref. 6 .) In the above equation, the 
summation is taken over the cyclic permutations, (i, j) — (1, 2), (2, 3), or (3, 1). 
Note that the displacement parameter 5 = xq cancels from this result also. 

In eq. (|28|l all parameters are determined from the two-particle correlation 
function with the exception of /3. This parameter is limited to the range of 
— 1 < (3 < 1, it is called the asymmetry parameter of Levy distributions, and can 
thus be determined from the relative momentum dependence of the normalized 
three-particle cumulant correlation function w. Note that symmetric stable 
distributions correspond to the case of (3 = 0, hence in that case </> = and 
w(l, 2, 3) = 1 even at large relative momenta. Thus all the essential parameters, 
a, f3 and 7 = R/2~ of stable source densities can be reconstructed from two and 
three particle correlation data. Let us emphasize that the result of eq. I)28|) is 
valid only within the plane-wave approximation and neglecting possible partial 
coherence in the source IIP) . 

In summary, we have determined the generic structure of Bose-Einstein cor- 
relations for the case when the coordinates of particle emission are obtained as 
convolution of many elementary subprocesses. Such processes can be attributed 
to parton-parton collisions, jet fragmentation, hadronization, rescattering and 
decay of hadronic resonances. 

Our choice of multiple convolution of various elementary probability laws was 
motivated by a recent paper by A. Bialas ^21, that considered Bose-Einstcin 
correlations for the case, when the radius of a Gaussian source fluctuates from 
event to event. Numerically, similar results were obtained by Utyuzh, Wilk and 
Wlodarczyk in ref. ^21 when considering Bose-Einstein correlations for sources 
with a fractal, power-law structure in space-time. 
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We find that the general shape of the Bose-Einstein correlation functions 
is a stretched exponential form. This is an exact analytic result, valid for all 
values of the relative momentum q, if the particle production is described by 
a stable law. Thus a new, experimentally measurable parameter is introduced 
to HBT or Bose-Einstein correlation studies, the Levy index of stability, a. 
The corresponding Levy stable distributions decay in the coordinate space as 
jjuj— i-a f or i ar g e values of \x\, see ref. [2] for greater details. 

Acknowledgments: T. Cs. would like to thank the Organizers of the 2nd 
Warsaw meeting on correlations and resonances for creating a meeting with an 
excellent and inspiring scientific atmosphere. This research was supported by an 
Alumni Initiatives Award of the Fulbright Foundation, an NWO - OTKA grant 
N25487, the Hungarian OTKA grants (T024094, T026435), T034269, T038406, 
T043514, the NATO PST.CLG.980086 grant and the MTA-OTKA-NSF grant 
INT0089462. 



References 

[1] R. Hanbury Brown and R. Q. Twiss, Nature 178, 1046 (1956). 



[2] T. Csorgo, S. Hegyi and W. A. Zajc, nucl-th/0310042 

[3] F. B. Yano and S. E. Koonin, Phys. Lett. B 78 (1978) 556 

[4] V. V. Uchaikin and V. M. Zolotarev, "Chance and Stability, Stable Distri- 
butions and Their Applications", VSP Science,1999, ISBN: 90-6764-301-7, 
596 pp. 

[5] V. M. Zolotarev, "One-dimensional Stable Distributions", Am. Math. Soc. 
Transl. of Math. Monographs, vol. 65, Providence, R.I. (Transl. of 1983 
Russian) 

[6] J. P. Nolan, Multivariate stable distributions: approximation, estimation..., 
http : // academic2 . american . edu/ ~ jpnolan/ stable/ overview . ps 

[7] J. P. Nolan, Fitting Data and Assessing Goodness-of-fit with Stable Distri- 
butions, Contribution to the Heavy Tails Conference, 
http : //academic2 . american. edu/ "jpnolan/ stable/DataAnalysis .ps 

[8] T. Csorgo, B. Lorstad, J. Zimanyi, Z. Phys. C71 (1996) 491 

[9] T. Csorgo and B. Lorstad, Phys. Rev. C 54, 1390 (1996) 

[10] T. Csorgo, |arXiv:hep-ph/0001233| Heavy Ion Phys. 15 (2002) 1-80. 

[11] T.Csorgo, B.Lorstad, J.Schmid-Sorensen, A.Ster, Eur. Phys. J. C9 (1999) 
275 

[12] A. Bialas, Acta Phys. Polon. B 23, 561 (1992). 

[13] O. V. Utyuzh, G. Wilk and Z. Wlodarczyk, Phys. Rev. D 61 (2000) 034007 



7 



